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Abstract: Piezoelectric ceramic Lead Zirconate Titanate (PZT) transducers, working on 
the principle of electromechanical impedance (EMI), are increasingly applied for structural 
health monitoring (SHM) in aerospace, civil and mechanical engineering. The PZT 
transducers are usually surface bonded to or embedded in a structure and subjected to 
actuation so as to interrogate the structure at the desired frequency range. The interrogation 
results in the electromechanical admittance (inverse of EMI) signatures which can be used 
to estimate the structural health or integrity according to the changes of the signatures. In 
the existing EMI method, the monitored structure is only excited by the PZT transducers 
for the interrogating of EMI signature, while the vibration of the structure caused by the 
external excitations other than the PZT actuation is not considered. However, many 
structures work under vibrations in practice. To monitor such structures, issues related to 
the effects of vibration on the EMI signature need to be addressed because these effects 
may lead to misinterpretation of the structural health. This paper develops an EMI model 
for beam structures, which takes into account the effect of beam vibration caused by the 
external excitations. An experimental study is carried out to verify the theoretical model. A 
lab size specimen with different external excitations is tested and the effect of vibration on 
EMI signature is discussed. 
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1. Introduction 

The EMI method has emerged as a widely recognized technique for dynamic identification and 
health monitoring of structural systems. The electromechanical admittance response of the smart 
system is derived from the dynamic interaction between the PZT transducer and the host structure. The 
EMI method has been proven direct and easy to implement. The method is typically applied using an 
electrical impedance analyzer, which scans a predetermined frequency range band in the order of tens 
to hundreds of kHz. In doing so, the complex admittance or impedance spectrum may be recorded and 
meaningful information containing the physical properties of the structure may be extracted. For SHM 
applications, these spectra can be compared at various times during the service lifespan of the 
structure, with which any change between the spectra is an indication of the presence of damage or 
material deterioration. 

Essentially the EMI method is founded upon the wealth of research accumulated from the modeling 
studies on electromechanical interaction between the structure and the PZT transducer, which started 
about two decades ago. Liang et al. [1] first proposed the EMI modeling approach for solving practical 
problems involving PZT actuator-structure interaction. This one-dimensional model has subsequently 
been developed and extended to two-dimension by other researchers [2-6]. The method was originally 
employed for the identification of structural parameters [7] and the design of piezoceramic driven 
electromechanical systems [8]. Sun et al. [9] first applied the EMI concept for the damage detection of 
a lab size truss structure, followed by Ayres et al. [10] who used the method for the detection of 
loosening prototype truss joints. Since then, the EMI based SHM method has been widely applied to 
various engineering structures, including aerospace [11, 12], civil [13-21] and mechanical [22, 23] 
structures. The EMI method has also been applied for vibration control [24]. 

However, the existing EMI method does not consider the structural vibration caused by the external 
excitation other than the PZT actuation, which means that the existing EMI method is only applicable 
to the static structures. Unfortunately, many structures are subjected to vibrations in practice. The 
vibrations of the structures will influence the results of SHM by causing the changes in PZT 
signatures. The changes in PZT signatures caused by the external excitation and the structural damage 
must be differentiated such that the correct information of structural health can be obtained. Therefore, 
this paper studies the effect of vibration on PZT signature by developing a new EMI model for beam 
structures with external excitations. Furthermore, experimental tests are carried out to verify the 
developed model. 

2. EMI Model 

As mentioned previously, the EMI method has been applied for the structures without excitations. 
However, no further research has been conducted for the EMI modeling of structures with excitations. 
This part will focus on the theoretical derivations of EMI models with excitations including both 
extensional and transverse vibrations. Uniform beams with simply supported boundary conditions will 
be considered. 

The equation of motion for the extensional vibration subjected to distributed force of an Euler- 
Bernoulli thin beam is given as: 
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E,W-p,AJi = f(x,t) 



(1) 



where E s is the elastic modulus; A s = h s b s is the cross-sectional area of the beam; p s is the material 
density of the beam; u{x,t) is the extensional displacement of the beam; f(x,t) is the harmonic 
distributed force; and b s and h s are the width and height of the beam, respectively. Assuming that both 

of the displacement and distributed force are variable separable: 



u(x,t) = U (x)e iM ; f(x,t) = f{x]e iwt ; and / (x) = ^F k cos—x 



(2) 



where F^, = — I /(jc)cos — ^J^j^is the frequency of external excitation; and l s is the length of the 

h 0 h 

beam. 

For a simply supported beam whose strains at two ends are zero, the boundary conditions are: 

u'(0,t) = 0 and u'(l s ,t) = 0 (3) 

From Eqs.(l) - (3), the extensional displacement for the forced vibration of the simply supported beam 
can be obtained as: 



i{x,t) = Y J - 



PA 



co 



=tCOS — xe 



Ps 



(4) 



The equation of motion for the transverse vibration of the same beam is given as: 

E,I,Y?+p,A,Y? = g(x,t) 



(5) 



where I s = h]b s l\2 is the moment of inertia of the beam; w(x, t) is the transverse displacement of the 
beam; the superscript iv means the fourth derivative with respect to x; andg(x,r)is the distributed 
force. Assuming that both of the displacement and distributed force are variable separable: 



w(x,t)=W(x)e i0} '; g(x,t) = g{x)e ic * and g(x)=JjG k sm™x 



(6) 



where G k = — | g (x) sin — xdx 

K 0 h 



For a simply supported beam whose displacements and bending moments at two ends are zero, the 
boundary conditions are: 

w{0,t) = 0;w"(0,t) = 0;w(l s ,t) = 0;andw"(l s ,t) = 0 (7) 

From Eqs. (5) - (7), the transverse displacement for the forced vibration of the simply supported beam 
can be obtained as: 



w(x,t) = Y J - 



M 



E S I S 2 
-co 



pA 



■ kn u 
sin — xe 

I 



(8) 
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Based on small displacement assumption, the following displacement field is defined for the 
extensional displacement of a generic point on the surface of the beam: 



h, 



( (x, t ) = u(x, t ) — - w'(x, t) 



(9) 



Therefore, the total extensional displacement is: 
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( kn^ 
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co 1 - 
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kn ,■„» khn ^ 

cos — xe' a — 

21. 
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k=l 



pA 



V h J av-\ 



EJ 2 

s s -co 2 



pAs 



k x io* (10) 



cos — xe 
I 



Hence, based on response compatibility, the total induced displacement of the PZT transducer is 
equal to the displacement difference between two discrete points on the surface of the beam, 
coincident to the ends of the transducer: 



u P2T - u s {x,t^ u s {x,t^ 
where x l and x 2 are the coordinates of the two ends of the PZT transducers, as shown in Figure 1 . 



(11) 



Figure 1. A simply supported beam bonded with a PZT transducer. 
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The following is the displacement of PZT vibration subjected to an electric potential: 



y ^M-^M) | K^F m (qt m ( ) -qf m {x 2 )) 



where the frequency of PZT actuation; </> n (x) = cos 



f \ 

fl7DC 



(12) 



(13) 



/ 



o„ = 



V *s J 



nn IE 
7 




( -J\ 
mm 



X=X2 



Sensors 2008, 8 



6850 



Q u = Q w = 0.5l s ; and F PZT is the magnitude of PZT actuation force. Detailed derivations of Eq.(13) can 
be found in [23]. 

Finally, we can obtain the electromechanical admittance of the PZT transducer: 

— sin(yl)ico PZT + 



Y = — = 



I Y*db 



h 



frsm 2 S{k) X {k)e^—h{co - co PZT ) 



i h 00 

m=l 



i(m-m PZT )t 



n=\ 

-diet) PZT I 



> + 



sim PZT lb 



(14) 



where: 



R 



y/ = ^cos(^Z) +— sin(^/) ; 



3{k) = 



khn 



g[^-Ol] 21, g[nl-^Y 



kn kn 
%(k) = cos— x 2 - cos— ^ ; 



S(n) = 
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Fjf is the complex Young's modulus; d is the piezoelectric strain coefficient; y is the wave number; 
R = K str K PZT ^ is the stiffness ratio; K str is the dynamic stiffness of the beam structure; 
K PZT = Y^bhr 1 is the quasi-static stiffness of the PZT transducer; s is the complex permittivity of 
PZT; and / , b and h are the length, width and height of the PZT transducer, respectively. 



3. Experimental Setup 

In order to verify the developed theoretical model, an experimental test is carried out on an 
aluminum beam bonded with one PZT transducer. The instruments and equipments used in this 
experiment are shown in Figures 2 and 3. 

The HP4192A Analyzer [25] is used to actuate the PZT transducer and simultaneously record its 
admittance signature. If multiple PZT sensors are used, a switch box is needed to make a multiple 
connections between the Analyzer and the PZT sensor. For this study, only one channel of the switch 
box is used since one PZT sensor is used in the experiment. In order to facilitate the automation of the 
testing, the program VEE Pro v. 6. 01 is used to control the Analyzer via a GPIB interface card installed 
in the PC. The scanning frequency range and resolution will be set with the Analyzer to record the 
PZT admittance signature. 
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Figure 3. External excitation system. 




The simply supported beam specimen is shown in Figure 4. The beam is stuck on the platforms by 
the blue tack. The mini-shaker generates external excitation to the beam at the location of 6cm from 
the left end of the beam. The PZT transducer is bonded to the beam at the central location. The 
properties of the PZT transducer and the beam specimen are listed in Tables 1 and 2, respectively. 
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Table 1. Properties of PZT transducer. 



Young's 
Modulus 

Y(GPa) 


Loss 
Factor 


Mass 
Density 

P (kg/m 3 ) 


Strain 
Constant d 

J 31 (m/volt) 


Permittivity 

£ 3 T 3 (farad/m) 


Dielectric 

Loss 
Factor 8 


Dimension 

(mm) 


Location 

(mm) 


68.9 


0.001 


7800 


-2.10E-10 


2.36E-08 


1.47E-02 


l=b=l0 
h=0.2 


100 mm 
from both 
ends 



Table 2. Properties of beam specimen. 



Length 

(mm) 


width 

(mm) 


Thickness 
(mm) 


Young's Modulus 
(GPa) 


Poisson's 
Ratio 


Mass 
Density 

(kg/m') 


Damping 
Ratio 


200 


30 


2 


68.9 


0.3 


2700 


0.005 



4. Results and Discussion 



As presented in the previous part, the theoretical derivation has been conducted to obtain the final 
result of the electromechanical admittance for a simply supported beam. And for practical cases, the 
transverse vibration of the beam is much more critical compared with the extensional vibration. 
Therefore, only the transverse vibration is considered in the experimental test. Then the 
electromechanical admittance of PZT can be simplified from Eq. (14) as: 

d 

d 
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-sm(yl)ico F 
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i(co — co pzl ) 



■ ldico c 



> + 



sico PZT lb 



(15) 



The above equation is applied to obtain the electromechanical admittance over the range of 
frequency bands. After substituting the values of various parameters into the theoretical formulation 
and a series of calculations by means of Matlab program, the theoretical results can be obtained. With 
the change in frequency and magnitude of the excitation force, the electromechanical admittance will 
change accordingly. Therefore, we will discuss the influence of these two factors. 

Figure 5 shows the theoretical prediction of PZT admittance for external excitations with the same 
magnitude but different frequencies, which are indicated at the bottom of the diagram. It is apparent 
that the conductance and susceptance of the PZT sensor are significantly affected by the excitation 
frequency. With the variations in frequency, both the peak location and the peak magnitude of the PZT 
admittance signature change. This result indicates that the frequency of external excitation is a critical 
parameter that influences the PZT signature. 

Figure 6 illustrates the theoretical prediction of PZT admittance for external excitations with the 
same frequency but different magnitudes. It can be observed that both the conductance and 
susceptance of the PZT sensor do not change drastically for different magnitudes of excitation force. 
Especially, close inspection shows that the peak locations of PZT admittance are almost the same for 
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different magnitudes of excitation. Moreover, with the increase of magnitude in excitation, the peak 
magnitude also increases. This result indicates that the magnitude of excitation force has significant 
impact on the peak magnitude but its effect on the peak location is negligible. 



Figure 5. Theoretical prediction of PZT admittance for different excitation frequencies. 
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Table 3. Theoretical and experimental peak locations for F=1.5sin (10007it) and F=1.6sin (1007tt). 



F= 


1.5 sin (10007rt) 


F= 


1.6 sin (lOOTrt) 




Pea 


c locations (kHz) 


Peak locations (kHz 




Theoretical 


Experimental 


Errors 


Theoretical 


Experimental 


Errors 


predictions 


results 




predictions 


results 




12.25 


11.80 


-0.45 


12.10 


11.95 


-0.15 


13.75 


14.25 


+0.50 


14.60 


14.60 


0 


18.50 


17.00 


-1.50 


15.90 


16.50 


+0.60 


21.30 


20.50 


-0.80 


21.15 


20.55 


-0.60 


24.75 


25.20 


+0.55 


25.60 


25.35 


-0.25 


27.25 


27.00 


-0.25 


26.70 


26.70 
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30.55 


30.15 


-0.40 


30.25 


30.30 


+0.05 


34.55 
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+0.20 
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-1.40 


34.85 
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-0.10 


38.50 


38.30 


-0.20 


36.90 


36.65 


-0.25 


40.00 


39.85 


-0.15 


38.25 


38.25 
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43.40 


42.50 


-0.90 


40.95 


40.00 


-0.95 


44.50 


43.75 


-0.75 


42.45 


43.65 


+1.20 


45.55 


45.00 


-0.55 


44.10 


44.10 
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46.70 
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45.00 


45.00 
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-1.25 
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To verify the developed EMI model, we first compared the predicted peak locations with the 
experimental ones. The peak locations for two different excitations are listed in Table 3. The 
differences between theoretical predictions and experimental results can be ascribed to several model 
assumptions as well as experimental errors. One of the major reasons for discrepancies could be that 
the effect of bonding layer (epoxy) between PZT and structure is not considered in the theoretical 
modeling. The errors between theoretical predictions and experimental results are also shown in Figure 
7. Overall, the accuracy of the theoretical model is acceptable. 



Figure 7. Errors between theoretical and experimental results. 
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Figure 8. Theoretical and experimental result of PZT admittance with external force 
F=1.5sin(10007it). 
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Figures 8 and 9 illustrate the theoretical and experimental admittance signatures. The comparisons 
show that the theoretical model is capable of predicting fairly accurately some of the 
electromechanical resonant peaks that can be found in the measured response. Nonetheless, 
discrepancies exist between the predicted response and the measured response. The most obvious 
discrepancy is found in the magnitude of the peaks. This discrepancy may be attributable to the 
simplified one-dimensional nature of the PZT transducer, as well as the assumptions made in the 
derivation of a closed-form solution for the dynamics of the structural substrate. Again, the accuracy of 
theoretical model is acceptable. 

Figure 9. Theoretical and experimental result of PZT admittance with external force 
F=1.6sin(100jrt). 
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In addition, we introduce damage in the specimen in the experimental test. This attempt is made to 
compare the peak shifts caused by damage and those caused by external excitation. The conductances 
for four cases (SI, S2, S3 and S4) are shown in Figure 10. SI and S2 represent the signatures of the 
specimen without damage while S3 and S4 represent those with damage. The damage in the specimen 
is induced by drilling a hole of diameter 5mm with distance of 4cm from the right end of the beam. 
There is no external excitation applied to the specimen for SI and S3 while the external excitation 
(F=sin (1000t)) is applied for S2 and S4. As expected, the shifts of the peaks are obvious (last 3 rows 
of Table 4). Comparison between SI and S2 shows the effect of external excitation. From Figure 10 
and Table 4, this effect is apparent and there is no consistent shifting trend in peaks. The peak shifts 
between SI and S3, and S2 and S4 represent the effect of damage on the signature. From Figure 10 
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and Table 4 (the last two rows), we can observe that the peak shifts due to damage are comparable for 
the specimen with or without external excitation. We can also see that the peak shifts due to damage 
and those due to external excitation are in the same order, implying that the effect of external 
excitation (vibration) on PZT admittance signature could not be neglected. For damage assessment, 
this effect should be eliminated; otherwise, it may lead to misinterpretation of the structural damage. 
Further study on how to eliminate/minimize the effect of vibration is definitely needed, which will be 
our future work. 




Table 4. Peak locations and shifts [unit: kHz]. 



Cases 


Peak locations and 


peak shifts 


SI 


11.5 


18 


22 


27.5 


30 


34 


35 


38.5 


41 


S2 


11.5 


17 


21.5 


27 


30 


33.5 


35.5 


39.5 


42.5 


S3 


12.5 


17 


21 


25.5 


30.5 


34.5 


36.3 


39.5 


42 


S4 


12.3 


16.2 


20.2 


25 


30.5 


34.2 


36.3 


40.2 


43.5 


S2-S1 


0 


-1 


-0.5 


-0.5 


0 


-0.5 


+0.5 


+1 


+1.5 


S3-S1 


+1 


-1 


-1 


-2 


+0.5 


+0.5 


+1.3 


+1 


+1 


S4-S2 


+0.8 


-0.8 


-1.3 


-2 


+0.5 


+0.7 


+0.8 


+0.7 


+1 



Figure 11 illustrates the effect of a broad band of frequency of external vibration on the 
conductance signature. The results are obtained using the EMI model developed in this paper. From 
Figure 11(a) where the external frequency <x> is much smaller than the PZT scanning frequency <x> PZT 
(5-50 kHz), co/co PZT < 1/25 , it can be observed that there is almost no peak shift in conductance 
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signature. Similarly, if co is much larger thanco PZT , e.g., col co PZT > 40 , the peak shift in conductance 
is also negligible, as shown in Figure 1 1(c). However, for 1/17 < col co PZT < 20 shown in Figure 1 1(b), 

obvious peak shifts can be observed. Together with other tries using the theoretical model, we 
conclude that if the external frequency is much smaller (e.g., co/co PZT < 1/25 ) or larger than 
(e.g., col co PZT > 30) the PZT frequency, the influence of external vibration on the EMI signature can be 
ignored. Otherwise, the effect of external vibration has to be considered in the EMI modeling. 



Figure 11. Theoretical prediction of PZT admittance for different excitation frequencies 
(F=0.5sinwt). 
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It is worth mentioning that we have focused on the peak location and peak shift because in the EMI 
technique, peak locations represent the natural frequencies of the host structure. Therefore peak shifts 
in admittance signature indicate the change of structural properties, which may be caused by damage 
or vibration. Other statistical method such as the root mean square deviation (RMSD) can also be used 
to analyze the signature [18]. 

4. Conclusions 

Although piezoelectric material has been successfully used for structural damage detection in many 
areas, the influence of the external excitations on the signature of PZT needs further study. This paper 
aims at ascertaining the effect of external excitation (vibration) on the PZT signature. The fundamental 
formulation of EMI model is first derived for a simply supported beam with external excitations. In the 
experiments, an aluminum beam specimen bonded with a PZT sensor is tested. The admittance 
signatures of the PZT sensor are measured for various excitations. The experimental results are 
compared with the theoretical predictions to investigate the reliability of the model. Overall, the 
theoretical model is capable of predicting the resonant peaks in the admittance signature. Further 
experimental tests are carried out to compare the effects of excitation and damage on the admittance 
signature. And the effective frequency range of external excitation is studied by the developed model. 
It is concluded that the effect of external excitation (vibration) on PZT admittance signature could not 
be neglected if the external excitation frequency and PZT's actuation frequency are comparable. For 
structural damage assessment, further study on how to eliminate this effect is needed. 
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